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Abstract. 

Fluctuations of conserved quantities in heavy-ion collisions have been argued to 
be diagnostic tools for the nature of the produced phase. These can be related 
to the predictions of quark number susceptibilities (QNS) from lattice QCD. Using 
the diagonal QNS, we extracted the Wroblewski parameter in a dynamical QCD 
computation. Our results on the cross correlations xbq, Xby, Xbs and xqy allow 
us to explore the charge and baryon number of objects that carry flavour. We present 
evidence that in the high temperature phase of QCD the different flavour quantum 
numbers are excited in linkages which are exactly the same as one expects from quarks. 
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1. Introduction 

As exciting experimental results keep pouring in from the Relativistic Heavy Ion Collider 
(RHIC) in BNL, New York, USA, the task of detecting the quark-gluon plasma (QGP) 
and establishing its various properties becomes more and more pressing. One needs 
therefore to look at as many signatures in as different a variety of aspects as possible. 
This has, of course, been going on both theoretically and experimentally. Fluctuations 
PP in conserved quantities such as the baryon number B, electric charge Q or strangeness 
S have been proposed as promising signals of QGP. Similarly, since very early days of 
relativistic heavy ion collisions, enhancement of strangeness production [2] has also been 
regarded useful indicator. 

The original proposals were usually based on simple model considerations as was 
the case for many other signals. The key fact employed there was that the up/down and 
strange quark masses are smaller than the expected transition temperature whereas the 
masses of the low lying hadrons, except for pion, are significantly larger. While these 
arguments on the differences between the two phases are qualitatively appealing, one 
has to face quantitative questions of details for any meaningful comparison with the 
data. Since the temperature of the plasma produced in RHIC, or even LHC, may not 
be sufficiently high for simple ideal gas pictures to be applicable, numerical estimates 
or indeed even the degree of utility of the proposal, could be misleading. Information 
obtained directly from the underlying theory, quantum chromodynamics (QCD) appears 
desirable, and the lattice formulation of QCD can potentially deliver that. 

A variety of aspects of the strangeness enhancement have been studied and many 
different variations have been proposed. One very useful way of looking for strangeness 
enhancement is the Wroblewski parameter [Hj. Defined as the ratio of newly created 
strange quarks to light quarks, 



the Wroblewski parameter has been estimated for many processes using a hadron gas 
fireball model [3]. An interesting finding from these analyses is that A s is around 0.2 in 
most processes, including proton-proton scattering, but is about a factor of two higher 
in heavy ion collisions. An obvious question one can ask is whether this rise by a factor 
of two can be attributed to the strangeness enhancement due to quark gluon plasma 
and if yes, whether this can be quantitatively demonstrated. We show below how quark 
number susceptibilities, obtained from simulations of lattice QCD, may be useful in 
answering these questions. 

An interesting property that the Wroblewski parameter possesses is that it is a 
robust quantity. Indeed, we [S] have recently argued that the ratios of susceptibilities, 
Ck/l, defined by 



2(ss) 



(1) 
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where \k and \l are quark number susceptibilities (QNS) for the conserved quantum 
numbers K and L and ok and <jl are the corresponding experimentally measured 
variances, are robust variables in the high T Phase. This is so both theoretically and 
experimentally, provided the two variances are obtained under identical experimental 
conditions and after removing counting (Poisson) fluctuations. We will show below our 
lattice QCD results in partial support of this claim. 

We employed such robust variables to address the important outstanding question 
of the nature of plasma (QGP) excitations. From several lattice QCD investigations, 
it has been known for a long time that the straightforward perturbation theory fails to 
describe the lattice data for equation of state for T > T c . Various resummation schemes, 
and phenomenological models have been tried to understand the nature of the QGP in 
this region. For T/T c > 3 — 5, lattice results on entropy density seem to be in agreement 
with such modified weak coupling pictures. Quark number susceptibilities provide an 
independent check on them, and yielded a strong support for these ideas. Screening 
mass determinations also suggested an agreement with an ideal Fermi gas of quarks for 
T > 2T C . Here we advocate a more direct approach. Creating an excitation of quantum 
number K, one can ask what else, e.g., like another quantum number L, does it carry. 
We address this by using the ratios of off-diagonal susceptibility xkl with \l '■ 

_ (KL)-(K)(L) 

C { KL)/L ~ {L2) _ {L)2 ■ (3) 

2. Quark Number Susceptibilities 

Let us first describe in this section the way one obtains various QNS from first principles 
using the lattice formulation of QCD. Assuming three flavours of quarks, and denoting 
by [if the corresponding chemical potentials, the QCD partition function is 

Z = jDU exp(-S G ) II BetM(m f ,fJi f ) . (4) 

f=u,d,s 

Note that the quark mass and the corresponding chemical potential enter only through 
the Dirac matrix M for each flavour. We employ the standard Wilson form for the 
gluonic action Sq and staggered fermions (with the usual square root method) to 
describe the quark matrix M. The formalism is, however, general and any change 
in either affects only the detailed expressions. The final physical results are expected to 
remain the same for sufficiently small lattice spacing a. 

Defining fi = fi u + \i d + fi s and /x 3 = \i u — /i^, the baryon and isospin densities and 
the corresponding susceptibilities can be obtained as: 

TdlnZ Td 2 \nZ 
Hl ~ V~cM' Xlj ~ ' (5) 

Similarly, Charge (Q), Hypercharge (Y), Strangeness (S) susceptibilities can be defined. 
QNS in (j3J) are crucial for many quark-gluon plasma signatures which are based on 
fluctuations in globally conserved quantities such as baryon number or electric charge. 
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Theoretically, they serve as an important independent check on the methods and/or 
models which aim to explain the large deviations of the lattice results for pressure 
P(/i=0) from the corresponding perturbative expansion. Here we will be concerned 
with i) extending our earlier quenched results on the Wroblewski parameter to full 
QCD with two dynamical quarks using our proposal |E] to estimate it from the QNS: 

A, = , (6) 

Xu "T Xd 

and ii) demonstrating the relevant degrees of freedom in QGP above T c . 

In order to use (JHJ) to obtain an estimate for comparison with experiments or for 
ii) above, one needs to compute the corresponding quark number susceptibilities on 
the lattice first and then take the continuum limit. All susceptibilities can be written 
as traces of products of M _1 and various derivatives of M with respect to fi. With 
m^i = m d , diagonal Xii s can be written as 

Xoo = ^[(0 2 (m u ) + iOnK))] (7) 
X33 = ^7 (0 2 (m u )) (8) 

X« = ^KP a (m.) + io al (mO>] • (9) 

Here 2 = Tr M~ X M'^ - Tr M" 1 M' u M~ l M' u , and O n {m u ) = (Tr M~ l M' u ) 2 , with 
the prime(s) denoting first(second) derivative of M(/Zj) with respect to /ij at /ij = 0. 
The traces are estimated by a stochastic method: Tr A = RiARi/2N v , and 

(Tr A) 2 = 2X^j=i(Tr A)i(Tr A)j/L(L — 1), where Ri is a complex vector from a 
set of N v , subdivided further in L independent sets. 

Higher order susceptibilities are defined by 
T d n \ogZ d n P 
f 9 '" V dfifdfXg- • - dfifdfig • - - ' 
These are Taylor coefficients of the pressure P in its expansion in /i. All of these can be 
written as traces of products of M~ l and various higher derivatives of M. These too can 
be evaluated using the Gaussian noise technique described above, but with increasingly 
larger number of vectors for progressively higher orders. One can build the expansion 
order by order to obtain information on the critical point of QCD in the T — fi plane, 
as done, for example in 

Figure 1 displays our results [5^ for the ratios Cb/q (left panel) and C(qy)/q (right 
panel) as a function of (square of) the lattice spacing a at 2T C , where T c is the transition 
temperature. The open symbols are for quenched QCD while the filled ones are for 
simulations with two flavours of light dynamical quarks. While the quenched results 
were obtained from a reanalysis of our earlier simulations [H], the data for the Nf = 2 
full QCD, also denoted as partially quenched to distinguish from the real world case 
of two light and one heavy strange quark, were obtained j7j on 4 x 16 3 lattices. The 
configurations were generated with a bare sea quark mass m = 0.1T C , which gives 
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Figure 1. Cb/q (left) and C(qy)/Q (right) exhibited as a function of l/N/~ oc a 2 
at T = 2T C , where a is the lattice spacing. Ratios of QNS are seen to be robust 
observables — being insensitive to both changes in lattice spacing a and the sea quark 
content of QCD. 



m w = 0.3m p . The setting of scale, the parameters employed and the statistics are 
detailed in Ref. [7j. To that set of data with T/T c = 0.75 ±0.02, 0.8 ±0.02, 0.85 ±0.01, 
0.9±0.01, 0.95±0.01, 1.00±0.01, 1.045±0.01, 1.25±0.02, 1.65±0.06 and 2.15±0.10, 
we added two more sets — 55 configurations separated by more than two autocorrelation 
times at T/T c = 0.975 ±0.010 (i.e., (3 = 5.2825) and 86 configurations, similarly spaced, 
at T/T c = 1.15 ± 0.01 (i.e., (3 = 5.325). 

In both the quenched and the partially quenched case and for both figures, the 
light valence quark mass m v , appearing in (JZJ)-©, is 0.03T C and the strange valence 
quark mass is T c . Both these ratios, and similar other ratios, display the property that 
r(a) = r ± 0(a n ) with n > 2, implying very good scaling behaviour. Furthermore, 
the results for the quenched and the partially quenched case on the N t = 4 lattice 
differ by a few per cent only, although T c in the two cases differ by a factor of 1.6. 
These two aspects of these ratios, very little dependence on the sea quark mass and the 
lattice spacing a, prompted us to term them as robust observables. Such robust ratios 
r can be reliably extracted on the smallest lattice already, which is very useful for the 
computationally expensive simulations with the dynamical quarks. One expects this 
robustness to persist at other temperatures above T c as well. At lower temperatures, all 
susceptibilities, including the off-diagonal ones, are of similar order and such robustness 
may not persist. It would, nevertheless, be interesting to test it carefully. 

3. Robust Predictions : Wroblewski Parameter 

Figure 2 exhibits some robust predictions of fluctuation measures from QCD, obtained 
from our dynamical simulations with light u and d quarks. Shown are the ratios Cx/q 
for the quantum number X indicated in the figure. The X = —S case is normalized by 
an extra factor of 2 to show it on the same plot. Note that the hierarchy seen in Figure 
2 is itself a robust feature. Our results indicate that experimental studies of Cs/q and 
Cb/q are most promising in terms of distinguishing between the two phases of QCD, 
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Figure 2. Some robust predictions of fluctuation measures from QCD. 



because they exhibit the largest changes in going from the hadronic phase to QGP. 

Recall that the strange quark susceptibility was obtained from the same simulations 
by simply choosing m v /T c = 1. Using (JHJ), A S (T) can then be easily obtained. Earlier 
we jHj had extrapolated A S (T) in quenched QCD to T c by employing simple ansatze. 
This became necessary in order to avoid any effects of the order of the phase transition; 
quenched QCD has first order transition while the Nf = 2 has a lot smoother transition. 
While the resultant A S (T C ) in quenched QCD displayed an impressive agreement [B] with 
the results obtained from the analysis of the RHIC and SPS data in the fireball model[3], 
clearly a cleaner determination with much better control of systematic errors needs its 
determination at T c in full QCD. Noting that A S (T) is a robust observable as well, one 
can do so from our N t = 4 simulations, at least for T > T c . 

Left panel of Figure 3 shows the Wroblewski parameter for full QCD as a function 
of T/T c . Its value at T c is seen to be X s ~ 0.4, in agreement with the value of extracted 
from experiments, when the freeze-out temperature is close to T c j3] . It is also a pleasant 
fact that at lower temperatures the ratio keeps decreasing. It is rather sensitive to the 
quark masses at those temperatures. The data sets A, D and C in the right panel of 
Figure 3 display xby / 'A^ s , i.e, the susceptibility \by divided by the square of the mass 
difference A us = m s — m u , for m s /T c = 0.1, 0.75 and 1 and m u /T c = 0.03. It is easy to 
show [Sj that such cross susceptibilities are zero in the limit of equal quark masses, and 
are oc A„ s . Strong kinematic effects are clearly visible near T c . It remains to be seen 
whether a realistic set of quark of quark masses, corresponding to the physical pion and 
Kaon masses, brings X s rapidly down by a factor of two or so below T c . 

The above mentioned nice agreement needs to be treated cautiously, however, in 
view of the various approximations made in equating the robust ratio determined on 
the lattice with the Wroblewski parameter extracted from the data. Let us list them in 
order of severity. 
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Figure 3. \ s as a function of temperature for QCD with two light dynamical quarks 
(Left panel) and various flavour symmetry breaking matrix elements (defined in the 
text) as a function of T/T c (Right panel). 



• The experiments at RHIC and SPS have nonzero albeit small /x whereas the 
above lattice result used /i = 0. Based on both lattice QCD and fireball model 
considerations, X s is expected to change very slowly for small \l. This can, and 
should, be checked by direct simulations. 

• Lattice simulations yield real quark number susceptibility whereas for particle 
production its imaginary counterpart is needed. Assuming that the characteristic 
time scale of plasma are far from the energy scales of strange or light quark 
production, one can relate [Hj the two to justify the use of lattice results in obtaining 
A s . Observation of spikes in photon production may falsify this assumption. 



4. Flavour Carriers : Quasi-quarks ? 

Unlike the gluons, quarks carry flavour such as electric charge or strangeness. Flavour 
in quark sector can potentially assist in identification of relevant degrees of freedom in 
QGP just above T c . One can look for correlations, by exciting one quantum number 
and looking for the presence of another. Choosing, baryon number and strangeness as 
these quantum numbers, Koch, Majumder and Randrup $| introduced the variable 

/-» or* 1 i Xus ~\~ Xds /-, i\ 

L'BS = —o^{BS)/S = 1 H ) l-L-U 

Xs 

in order to distinguish between bound state QCD -sQGP [ID] - and the usual picture 
of the excitations in the plasma phase of QCD. This is expected to have a value of 
unity if strangeness is carried by (ideal) quarks, since 5 = 1 always comes linked with 
B = —1/3. In [9^ it was shown that bound state QGP gives a value of Cbs ~ 2/3 (for 
T > T c ). Charge and strangeness correlation offers another similar possibility of being 
unity, if strangeness is carried by quarks : 

n Qri 1 ^Xus — Xds / 1n \ 
L,Q S — 0O(Q5)/ 5 - 1 . [LZ) 
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Figure 4. Cbs and Cqs, as functions of T/T c . The quark masses used are 
ni-ud = 0.1T C and m s = T c , although in the high temperature phase there is no 
statistically significant dependence on the quark masses. 



In Figure 4, we show the first lattice QCD results 5j on Cbs an d Cqs from our 
dynamical two light quark simulations. Note that while both are different from unity 
below T c , they become close to unity immediately above T c , strongly indicating that 
unit strangeness is carried in QGP by objects with baryon number -1/3 and charge 1/3 
very close to T c and beyond. The trends for T < T c as well as the order of magnitude 
estimated in the hadron gas models are in agreement with our lattice results. We 
varied the strange quark mass, m s /T c , between 0.1 and 1.0 and found that for T > T c it 
does not alter the value, ~ 1, or the visible T-independence. A natural explanation of 
the T-behaviour is provided if strange excitations with nonzero baryon number become 
lighter at T c . Furthermore, the T-independence suggests dominance of a single such 
excitation. Similar results were also obtained in j3] in the light quark sector, from e.g., 
C(bu)/u an d C(Quyu, linking w-flavour carriers to B = 1/3 and Q = 2/3 objects. 

A possible interpretation of the observed very small deviations from unity of Cbs 
and Cqs is as follows. Colour interactions merely dress up quarks. These effects can 
be computed in weak coupling limit, i.e, at high T, where the flavour linkages of quarks 
are natural. Close to T c , the flavour linkages remarkably continue to persist as before, 
but they are no longer computable as the coupling is presumably not weak. This leads 
us to term them therefore as quasi-quarks. 

5. Summary 

Quark number susceptibilities which can be obtained from first principles using lattice 
QCD contribute substantially to the physics of RHIC signals. Since the ratios of quark 
number susceptibilities, Ca/b were shown to depend rather weakly on the lattice spacing 
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and the sea quark content of QCD in the high temperature phase, they are theoretically 
robust variables. Consequently, one can compute them reliably relatively easily on small 
lattices, as we did. Being ratios, experimentally they are likely to be free of various 
systematic errors and thus again robust experimental observables. We presented first 
full (two light dynamical quarks) QCD results for the Wroblewski Parameter \ S (T). 
Near T c , these are found to be in agreement with the RHIC and SPS results. Being a 
robust observable, small lattice cut-off effects expected in A s . 

We demonstrated that the high temperature phase of QCD essentially consists of 
quasi-quarks by exploiting the flavour linkages of quarks. In particular, we showed that 
unit strangeness is carried by an object with baryon number —1/3 and charge 1/3, as 
seen in Figure 4. Moreover, this correlation does not depend on the strange quark mass 
even when it is as large as T c . Similarly, in the light quark sector one finds that u and 
d quantum numbers are not produced together, and that the u flavour is carried by 
excitations with baryon number +1/3 and charge +2/3, whereas the d flavour is carried 
by particles with baryon number +1/3 and charge —1/3. 
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